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ABSTRACT 

Wc study the visibility of the Lya emission line during the Epoch of Reionization 
(EoR). Combining galactic outflow models with large-scale semi-numeric simulations 
of reionization, we quantify the probability distribution function (PDF) of the fraction 
of Lya photons transmitted through the intergalactic medium (IGM) , Tigm ■ Our study 
focusses on galaxies populating dark matter halos with masses of Mhaio = 10^°Mq at 
z = 8.6, which is inspired by the recent reported discovery of a galaxy at z = 8.6 with 
strong Lya line emission. For reasonable assumptions, we find that the combination 
of winds and reionization morphology results in Tigm >10% [50%], for the majority 
of galaxies, even when the Universe is ^ 80% [60%] neutral by volume. Thus, the 
observed strong Lya emission from the reported z = 8.6 galaxy is consistent with a 
highly neutral IGM, and cannot be used to place statistically significant constraints 
on the volume averaged neutral fraction of hydrogen in the IGM. We also investi- 
gate the implications of the recent tentative evidence for a observed decrease in the 
'Lyman Alpha Emitter fraction' among drop-out galaxies between z = 6 and z = 7. 
If confirmed, wc show that a rapid evolution in xhi will be required to explain this 
observation via the effects of reionization. 

Key words: galaxies: high redshift - galaxies: stellar content; cosmology: dark ages, 
reionization, first stars - cosmology; early Universe - cosmology; diffuse radiation - 
cosmology; large-scale structure of Universe - radiative transfer 



1 INTRODUCTION 

The Wide Field Camera 3 on board of the Hubble Space 
Telescope has enhanced our ability to observe galaxies at red- 
shifts great than six, so- far obtaining ~ 100 likely candidate 
galax i es at z =7, 8 (e.g. 



I2OIOI : iFinkelstein et all 



Bouwens et al.ll2010al: [Bunker et al.l 



2OIOI : lYan et all bOld ). The James 



Webb Space Telescope is expected to probe galaxies a few 
magnitudes deeper, and also to spectroscopically confirm 
the redshifts of the existing candidates. One of the key pre- 
dicted properties of young, metal poor galaxies in the high- 
redshift Universe are prominent nebular emission lines, dom- 
inated by hydrogen Lya (A = 1216 A, see e.g. Johnson et 
al. 2009, Pawlik et al. 2010). The first generation of galaxies 
are likely to have been strong Lya emitters, wi th equiv- 
alent widths possibly as high as EW~ 1500 A (|Schaererl 
I2002l . l2003l : I Johnson et al I2OO9I . also see Partridge & Peebles 
1967). 



During the epoch of reionization (EoR), the Lya emis- 
sion line may be difficult to observe, due to the large opacity 
of the intervening neutral intergalactic medium: for exam- 
ple, a source needs to be embedded in a ^1 Mpc HIT re- 
gion to allow Lya photons to redshift far away from the 
line center before they reac h the IGM (e.g.. iMiralda-Escudd 
ll998l : [Cen fc Haimanll2000l '). For the galaxy to generate such 
a large HII region, its ionizing luminosity would have to be 
unphysically large (unless there is a bright quasar associated 
with the galaxy, see Cen & Haiman 2000). However, several 
effects have been shown to boost the detectability of the Lya 
flux: ( i) source clust ering, which boosts the sizes of HII re- 
gions (iFurlanetto et al. 2004; Mcs inger fc Furlanettoll2008bl : 
iMcQuinn et all 120071 : ilhev et all boOSl ) iii) the patchiness 
of reionization, which may give rise to significant fluctua- 
tions in the IGM opacity between dif ferent sightlines, as well 
as a steeper absorption prof ile (e.g.. iMesinger fc Furlanettd 
l2008al :[M cQuinn et al.ll200a ). and {iit) radiative transfer ef- 
fects through outflows of interstellar (ISM) H I gas, which 
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can impart a redshi ffl to the Lyg photons before they 
emerge from galaxies (|Santosl2004l :l iDiikstra fc Wvithel2O10l . 



but also see Barne s et al . 2011). 

iLehnert et ai] (|201(]| ') recently reported a detection of 
strong Lya line emission from a Y105 drop-out galaxy in 
Wide Field Camera 3 obseryations of the Hubble Ultra 
Deep Field. The Lya line implies that the galaxy is at 2 ~ 
8.56, which is the highest redshift of any spectroscopically- 
confirmed object to date. Interestingly, when taken at face 
value, the observed Lya line is strong, with an observed 
equivalent width (EW) of ~ 100 A (see §111). Motivated by 
this observation, we study the visibility of the Lya emission 
line during the EoR, and compute the total fraction of emit- 
ted Lya photons that the IGM transmits directly to the ob- 
server, TiGM- We simultaneously include the inhomogeneous 
large-scale reionization morphology, peculiar velocity offsets 
of the galaxies and IGM, and radiative transfer through the 
galactic outflows that is calibrated by ob servations of Ly- 
man Alpha Emitters (LAEs) a t z < 6 (|Verhamme et all 
l2006l . l2008l : IVanzella et al]|2010h . This, in combination with 
the fact that we compute the full 7igm-PDF, clearly distin- 
guishes our analysis from previous work. 

Our models do not include dust, which at the red- 
shifts of interes t (z = 7 — 9) is likely a good a p proximation 
(•Stanwav ct al.' '2005'; iBouwens etaP boiObI : faaves et all 
[2011.; .Blanc ct al. .201l1). The dust opacity to Lya photons 
inside galaxies is expected to increase towards lower redshift, 
as the cumulative dust content of the Universe increased 
with cosmic time. This expected evolution has a different 
sign than the IGM opacity which decreases towards lower 
redshift. Thus, our discussion regarding the redshift evolu- 
tion of the IGM opacity is likely conservative. 

The outline of this paper is as follows: in §[2] we describe 
our models for galactic outflows (i )2.1|l and IGM opacity 
f i]2.2|l . In §[3]we present the corresponding Lya transmission 
fractions. Within this context, we interpret the observations 
of Lya emitting galaxies at 2 > 6 in § |4] We compare our 
results with previous work in § (5) Finally, we present our 
conclusions in § |6l We adopt the background cosmological 
parameters (Oa, ^m, i^b, n, as, Ho) = (0.72, 0.28, 0.046, 
0.96, 0.82, 70 km s"^ Mpc~^ ), matching the five- year re- 
suhs of the WMAP satellite (|Komatsu et al.1 120091 '). Unless 
stated otherwise, we quote all quantities in comoving units. 



2 THE MODEL 

As mentioned above, our model has two components; (i) the 
intrinsic Lya line, which has been processed by Lya scat- 
tering through galactic outflows, and (ii) the IGM opacity, 
which further processes the line through both resonant and 
damping wing absorption. We describe these below in turn. 



^ The peculiar velocity of a galaxy can also redshift Lya photons 
away from resonance before they escape into the surrounding in- 
tergalactic medium (Cen et al. 2005). However, this redshift is 
typically significantly smaller than the redshift imparted by galac- 
tic outflows (see Dijkstra & Wyithe 2010). 



i20 



200 km s" 



lies in the middle 



2.1 Modeling Galactic Outflows 

We employ the wind models described in Dijkstra & Wyithe 
(2010) . In these models the outflow is represented by a spher- 
ically symmetric shell of HI gas that has a column density 
A'hi, and outflow velocity t'wind- The shell has a thickness 
of 0.1 kpc and a radius of 1.0 kpc (proper, also see Ahn 
et al. 2003). The outflow surrounds a central Lya source 
which emits photons at line center. Since we focus solely 
on the total fraction of Lya photons transmitted to the ob- 
server, the total Lya luminosity of the source is irrelev ant. 
A Monte-Carlo code (described in IDiikstra et al.ll2006al ) ac- 
curately follows the propagation of Lya photons through 
the optically thick outflow. These models are very simi- 
lar to the models that reproduce observed Lya line pro- 
files at 2 = 3 - 5 (|Verhamme et al.l |2006| , |2008| . see Dijk- 
stra & Wyithe 2010 for a discussion of caveats etc.). In this 
paper, we focus on wind models with A'hi = 10^*^ cm~^ 
and A^Hi ~ 10^^, and wind velocities of w„ind = 25 km s~^ 
and Uwind ~ 200 km s~^. For comparison, Verhamme et 
al (2008) reproduced observed Lya line shapes for LAEs 
at 2: < 5 with 25 km/s^ fwind ^ 400 km/s, and that 
2 X 10^^ cm" 
(iVni = 10^" cm-^ t;„ind 
of this range. 

The impact of galactic winds on the Lya radiation field 
depends on the covering factor of the galactic outflow: for 
low covering factors only a small fraction of the emitted 
Lya photons will be Doppler boosted to frequencies where 
the IGM opacity is reduced. The issue of the wind covering 
factor is therefore related to the overall impact of galac- 
tic winds on the Lya radiation field. Existing observations 
(which extend out to 2 ~ 6) indicate that winds play an im- 
portant role in the scattering process, implying a large cov- 
ering factor of the outflowing scattering material (we discuss 
this further in §[5]). Of course, it remains a possibility that 
winds properties, such as their covering factor, were differ- 
ent at 2 > 6. This would affect the impact of the (especially 
ionized) IGM on Lya emission lines, and hence the observed 
redshift evolution of LAEs at these redshifts (see i; 14. 2p . 



2.2 Modeling the IGM Opacity 

We use the publicly-available, semi-numerical code DexlvQ 
to generate evolved density, velocity, halo, and ionization 
fields at 2 = 8.56. This code and detailed tests are presented 
in Mesinger & Furlanetto (2007), Mesinger et al. (2010) and 
Zahn et al. (2010), to which we refer the reader for details. 
Here we briefly summarize our simulation. 

Our simulation box is L = 250 Mpc on a side, with 
the final density, peculiar velocity, and ionization fields hav- 
ing grid cell sizes of 0.56 Mpc. Halos are filtered out of 
the 1800'' linear density field using the excursion-set formal- 
ism, and then mapped to Eulerian coordinates at 2 = 8.56 
through perturbation theory (Zeldovich 1970) . Perturbation 
theory is also used to generate the evolved density and pe- 
culiar velocity fields. Corresponding ionization fields are cre- 
ated according to the excursion-set prescription described in 
Mesinger & Furlanetto (2007), with the modification from 



http: / /www. astro. princeton.edu/~mesinger/Sim. html 
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Zahn et al. (2010) to account for partially ionized cells. This 
prescription compares the number of ionizing photons pro- 
duced in a region of a given scale to the number of neu- 
tral hydrogen atoms inside that region. We generate a suite 
of ionization fields at various values of xhi by varying the 
ionization efficiency of sources assumed to be hosted by 
atomically-cooled halos, with masses of Mhaio > 10* M0. All 
of these fields have been extensively tested against hydrody- 
namical cosmological simulations, and good agreement was 
found well past the linear regime (Mesinger & Furlanetto 
2007; Zahn et al. 2010; Mesinger et al. 2010). 

We then extract ~ 10* line-of-sights (LOSs) centered 
on halos in the mass range, IO^^Mq < Mhaio < 3 x 
10^° Mq. This choice of host halo masses is motivated by 
the UV derived star formation rate (SFR) of 2-4 Mq yr"^ 
(jLehnert et al.ll2010l l. which corresponds to a halo mass of 
~ lO^^M© in the cosmological hydrodynamic simulations of 
Trac & Cen (20070. Opacities at wavelengths surrounding 
the Lya line are computed for each LOS, integrating under- 
neath a Voigt absorption profile (e.g. Rybicki & Lightman 
1979), and including contribution from both the ionized and 
neutral IGM, out to distances of 200 Mpc away from the 
source. Inside HIT regions, a residual HI fraction is computed 
assuming ionization equilibrium with an ionization rate of 
Thii = 0.5 X 10~^^ s~^, in rough agreement with estimates 
obtained from the z -^5-6 Lya forest (e.g.. Fan et al. 2006; 
Bolton & Haehnelt 2007j3. We include the peculiar veloc- 
ities of both the source ha lo and the absorbing gas, which 
can be very important (e.g . iDiikstra et al.ll2007bl : lniev et al] 
l2008l : lLaursen et al.|[201lh . 



3 RESULTS 

To quantify the detectability of the Lya line, we compute 
for each LOS the total fraction of emitted Lya photons that 



^ Our results are not sensitive to uncertainties of a factor of few 
in the host halo mass, because the halo bias evolves relatively 
slowly with mass in this range (Mesinger & Furlanetto 2008a; 
McQuinn et al. 2008; Mesinger & Furlanetto 2008b). 
* We emphasize that there is currently no evidence that reioniza- 
tion has completed at these redshifts (Lidz et al. 2007; Mesinger 
2010; McGreer et al. (2011). If the observed quasar spectra go 
through regions of pre-ovcrlap neutral gas, the inferred value of a 
homogeneous Ftot would include contributions from the neutral 
IGM (F ~ 0). Therefore, the derived values of Ftot (from, e.g.. 
Fan et al. 2006; Bolton & Haehnelt 2007) can be treated as lower 
limits for the ionization rate inside the ionized component of the 
IGM, Fhii. 
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Figure 1. This Figure shows an example of a Lya lino profile. 
The solid line shows the Lyo spectrum-in units of the peak flux 
dcnsity-that emerges from the galaxy after the photons have scat- 
tered through the HI outflow. The blue dotted line shows the IGM 
opacity, exp[— -riGM('^)] as a function of velocity offset Av for 
5;hi = 0.51. The red dashed line shows the spectrum of photons af- 
ter processing the flux through the IGM, i.e. Jiv) exp[— riQM('^)]- 
This Figure nicely shows the impact of resonant scattering in the 
ionized IGM at Av ^170 km s~^, and the impact of the damping 
wing optical depth at Av ^170 km s~^. The latter varies only 
weakly with frequency. 



the IGM transmits directly to the observer, Tigm, a^ 

TiGM = / di'J{iy)exp[-TiGM{v)], (1) 



where tigm(j^) is the optical depth of the intervening IGM at 
frequency v (computed as described in H2.2p . and J(i') is the 
normalized (i.e. J J{v)di' = 1) Lya spectrum of Lya pho- 
tons that emerges from the galaxy (computed as described 
in H2.ip . The integral runs from ~ 10^ km s~^ blueward to 
~ 3 X 10'' km s~^ redward of the Lya line resonance, which 
spans the full range of velocities that is covered by the Lya 
profile that emerges from the galaxy. The quantity Tigm is 
also referred to as the 'IGM transmission fraction'. 



3.1 An Example Line Profile 

An example on an observed Lya line profile is shown in Fig- 
ure [T] In this plot, the solid line shows the Lya spectrum-in 
units of the peak flux density-that emerges from the galaxy 
after the photons have scattered through the HI outflow (in 
this case Wwind = 200 km s~^, and A'hi = 10'^'' cm~^). 
Most of the line flux is systematically redshifted relative 
to the galaxy's systemic velocity. The flux density peaks 



^ Photons that are scattered in the neutral IGM produce diffus e 
Lya halos around individual sources (e.g. IZheng et al]|2010d) . 
This emission is several orders of magnitude fainter than the de- 
tection threshold of the deepest observations to date (see Dijk- 
stra & Wyithe 2010). Lya radiation that is resonantly scattered 
in close proximity ( <;10 kpc) to the ga laxy can give rise to a 
brighter Lya halos ll Zheng et al.l l2010alf^. However, winds re- 
duce the brightness of these halos llDiikstra fc Wvithell2oTol) . In 
any case, by ignoring this resonantly scattered component, we un- 
derestimate the true Lya flux that can be detected from galaxies 
during the EoR, which only renders our results conservative. 
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Figure 2. Left panel: This Figure shows the probabihty density, dTo^^^J^' that the IGM transmits to the observer a fraction of emitted 
Lya photons in the range logTioM i dlog7iGM/2, for galaxies populating dark matter halos of 10^'' Mq < Mhaio < 3 X 10^'' Mq in 
a universe with a neutral fraction of ijji = 0.91 (by volume) at z = 8.6. The red dashed line (black solid line) shows the model with 
''wind = 25 (200) km s~^. In both models A^hi = lO'^" cm"'^. The blue dotted line shows the 'no-RT' model (see text). This Figure 
illustrates that (i) the IGM can transmit a significant fraction of Lyo photons, despite the fact that reionization has only just started 
(text), and (ra) the IGM becomes even more 'transparent' when winds are affecting Lyo scattering in the ISM. Right panel: Same as the 
left panel, but now we plot the cumulative distribution function (CDF), P{> Tigm)- We find for example that Tiqm > 10% for ~ 10% 
of all halos in the 'no-RT' model, and that this fraction is boosted when winds are present. 
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Figure 3. Same as Figure[2] but with an enhanced HI column density in the wind of A'hi = lO'^^ cm~^. Frequency diffusion that occurs 
as Lyo photons scatter through an extremely opaque medium causes a larger fraction of photons to emerge at larger redshifts from the 
line line center, which reduces the effective IGM opacity. 



at ~ 2t;wind, which is expected for radiation that scattered 
to the observer from the outflowing gas on the far side of 



to tne ODserver tror n the outnowmg gas on the tar side oi 
the Lyg source (seelAhri et al.ll200a : IVerhamme eraLll2006l : 
iDiikstra fc Wvithell2010l )~ 

The blue dotted line shows the IGM opacity, or more 
precisely exp — [tigm{i^)] as a function of velocity off-set 
Aw for xm = 0.51. In this particular example infalling 
(ionized) gas provides a large opacity to Lya photons, 



frame of the galaxy, e.g. 


Santosll2004l: iDiikstra et al.ll2006bl. 


2007bl; Ilicv ot al.' '2008 


; 'Daval et al.l 201ll: Zheng et all 


2010al:lLaurscn ct al. 2011). That is. ttgm » 1 at Au <-I-170 



km s~^. At redder wavelengths there is no gas that falls in 
fast enough for the Lya photons to appear at resonance. 
At these frequencies {Av Jjl70 km s~^) the IGM opacity is 
dominated by the damping wing optical depth of the neutral 
IGM which is a smooth function of frequency. 

The red dashed line shows the spectrum of pho- 
tons after processing the flux through the IGM, i.e. 
J(i/) exp[— tigm(i^)]. For this particular example, a fraction 



Tigm ~ 0.57 of all photons is transmitted directly to the 
observer. For comparison, had we assumed that all photons 
emerged from the galaxy with a Gaussian emission line, cen- 
tered on the galaxy's systemic velocity and with a standard 
deviation of cr = Vdrc ~ 80 km s~^, we would have found 
that Tigm = 0.01. 



3.2 The Probability Distribution Function (PDF) 

of Tigm 

In Figure[2l we show the PDF, diog^j ^ {left panel) , and the 
cumulative distribution function (CDF), P(> Tigm) {right 
panel), for xm ~ 0.91. The red dashed lines {black solid 
lines) correspond to models with i^wind = 25 (200) km s~^, 
and JVhi = 10^° cm"^ (both models). For the low (high) 
wind velocity we find that the log Tigm-PDF peaks around 
~ 0.06 (0.1). Increasing the wind velocity clearly shifts the 
Tigm-PDF to larger values. As the Lya photons scatter off 
of the receeding outflows, they are Doppler shifted to larger 
effective redshifts (e.g., Ahn et al. 2003, Verhamme et al. 



© 2009 RAS, MNRAS OOO.nifTOl 



The Detectability of LAEs during the EoR 5 




Figure 4. Same as Figure [2] but witli a reduced global neutral hydrogen fraction, xhi = 0.80. The shifts of the Tigm-PDFs are due to 
the reduced neutral hydrogen content of the Universe. 
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Figure 5. Left panel: The fraction of hales with 10 Mq < Mhalo < 3 X 10 A/q for which 7igm > 0.1 as a function of 2'hi, assuming 
''wind = 25 (200) km s denoted by red filled circles (blue filled squares). In both models Afjji = 10^" cm"'^. Right panel: same as the 
left panel, but now we plot the fraction of halos for which Tigm > 0.3 (top set of points), and Tigm > 0.5 (bottom set of points). This 
panel shows that the fraction of halos with Tigm > 0.3 only becomes less that ~ 10% when xhi ^0.90. We therefore conclude that the 
observed strength of the Lya line in the z = 8.56 galaxy is not surprising, unless the Universe was more than 90% neutral. 



2006). Therefore, larger wind velocities cause a larger frac- 
tion of Lya photons to emerge at frequencies where they 
are immune to the opacity in the ionized infalling IGM (see 
Fig. [T]). Furthermore, this Doppler-shifting from the winds 
means that by the time the Lya photons reach a neutral 
patch of the IGM, their absorption cross-sections are fur- 
ther out on the damping wing tail. Therefore, outflows re- 
duce the impact of both resonant and damping with opac- 
ities (see Fig [1]). Specifically, the right panel shows that 
TiGM > 10% for ~ 30% (~ 60%) of LAEs for Uwind = 25 
km s~^ (I'wiiid = 200 km s~^). 

To underline the effect of winds, we compare to a model 
in which we only evaluate the damping wing optical depth, 
td at line center, i.e. we set Tigm = — lnrD(Az/ = 0). This 
model is referred to as the 'no-RT' model, as it corresponds 
to a model in which no scattering of Lya? photons occurs 
in the either the ISM or the ionized IGM, and is repre- 
sented by the blue dotted line. The Tigm is clearly skewed 
more to lower values for this no-RT ('RT' stands for ra- 
diative transfer) model. That is, without winds a neutral 
IGM dramatically attenuates the transmission of the Lya 
Une. Only ~10% of LAE have transmission fractions greater 
than 0.1 (as has been demonstrated repeatedly in the past, 
e.g.. Gen & Haiman 2000, Santos 2004; Furlanetto et al. 



2004; McQuinn et al. 2008; Mesinger & Furlanetto 2008a,b). 
Winds therefore clearly boost the detectability of Lya emis- 
sion from galaxies during the early phases in the EoR. 

Figure [3] shows the same quantities as Figure [2l but 
for wind models with A^'hi = 10^^ cm"^ The 'no-RT' 
model is of course unchanged. The Tigm-PDFs are shifted 
to larger values, because resonant scattering through very 
opaque media (the line-center optical depth to Lya photons 
is TO = 5.9 X 10^(AfHi/102' cm-2)(rgas/10'' K)-^''^), results 
in frequency diffusion which increases with optical depth tq 
(e.g. Harrington 1973, Neufeld 1990). As a result of this fre- 
quency diffusion, a larger fraction of photons will emerge 
at larger redshifts from the line line center, which further 
reduces the effective IGM opacity. 

Figure |4] shows the same as Figure [2] but for a lower 
volume averaged neutral fraction of xm — 0.80. We find 
that Tigm > 10% for ~ 35% for the no-RT model, and that 
TiGM > 10% for ~ 75% (~ 95%) of all halos for the wind 
model with Wwind = 25 km s~^ (t'wind = 200 km s~^). This 
shift of the Ti:GM-PE)F arises because of the reduced neutral 
hydrogen content of the Universe. 
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4 COMPARISON TO RECENT DATA 

4.1 Interpretation of the Recent Observations of a 
2 = 8.56 Galaxy 

The observed Lya luminosity of tlie z = 8.56 galaxy is 
La = [5.5 ± 1.0 ± 1.8] X 10"^ erg s'^ (the first number de- 
notes the 1 — cr-uncertainty, while the second denotes the 
systematic uncertainty), while the UV luminosity density is 
L^(A = 1700 A) = io28 3±o-2 erg s"^ Hz^^ at A = 1700 
A (rest-frame, Lehnert et al. 2010). From these observed 
strengths of the line and continuum, we can constrain the 
observed Lyg rest fram e equivalent width (REW) to be 
jPiikstra fc Westra|[201ol ): 



REW 



Lx{X — \a) Vol Li,[\ — Aq) 

'1216\/3-2 



(2) 



Vo L^(A = 1700 A) 



/1216y 
5) V 1700/ 



= 1361??, ± 45 A. 



We assumed that the UV continuum slope is /3 = 2, which is 
appropriate for s tar forming galaxi e s with strong Lya emis - 
sion aXz = 3-7 (jStark et al.ll2010al '). lBouwens et al](j2010bl ) 
determined /3 ~ 3 for candidate 2 ~ 7 (zsso drop-out) galax- 
ies, with Muv ~ — 19 to ~ —18. Inserting = 3 into Eq 2 
reduces the expectation value for REW to REW= 97 A. 

The quoted 1 — a uncertainty on REW is dominated by 
the uncertainty in the observed UV flux density. The uncer- 
tainty on the observed REW is large, but the expectation 
value is remarkably large (even for /3 = 3). The intrinsic 
Lya REW, REWint, depends quite strongly on gas metal- 
licity, and whether the galax y is forrning s tars in a burst 



Y ,, , 

or continously (e.g. Fig 7 of ISch acrcr' 2003), and whether 

the 'case-B' approximation is valid (Raitcr ct al. 2010l). The 
m aximum possible va lue appears to be REWmax = 3000 
A (|Raiter et all I2OIOI '). This maximum value is reached if 
this galaxy formed stars from metal-free gas in a burst 
with a top-heavy initial mass function. Under this assump- 
tion, the data requires that the IGM transmits Tigm ~ 



0.04 



-F0.03 . 
-0.02 • 



/ REWi„t \ ' 
V 3000 A / 



, where REWint denotes the intrinsic 



(or emitted) REW. In theory, we could explain the observed 
REW for TiGM >0.04tg °^. However, throughout we focus 
on the more conservative requirement Tigm > 0.1, which 
requires smaller values for the intrinsic REW. 

For comparison, Figures OH show that even with low 
velocity outflows (Wwind = 25 km s"^) TiGM > 0.1 for a 
significant fraction LAEs, even if the Universe is ~ 91% 
neutral. This is illustrated more explicitly in Figure[51 where 
the left panel shows the fraction of LAEs for which Tigm > 
0.1 as a function of xm for Uwind = 25 (200) km s ~^ as red 
filled circles {blue filled squares), where A*'hi = 10^" cm~^ in 
both models. This fraction is ^50% for both models when 
xhi ;^80%. The right panel is the same as the left panel, 
but for Tigm > 0.3 (upper set of points), and Tigm > 0.5 
(lower set of points). This panel shows that the fraction of 
LAEs with Tigm > 0.3 only becomes less than ~ 10% when 
xm ^0.80. We therefore conclude that the strength of the 
Lya line of the = 8.56 galaxy is not surprising, unless the 
Universe were more than 90% neutral by volume. 

This latter statement is quantified in Figure [6l where 
we plot the conditional probability distribution for xm, 
if we require that Tigm > 0.1 (the three lines rep- 
resent the same models that were shown in Fig 2-4). 
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Figure 6. This figure shows the marginahzed PDF for xm, given 
that TiGM > 0.1, denoted by ''^^^"'^^^f ■ Different lines 

represent the different models that were also shown in Figures 
2-4. This Figure shows clearly that it is difficult to rule out even 
large values of xm with confidence— especially when winds are 
included in the modeling (see text). 

This PDF follows from Bayes theorem afl ^(ShiITigm) = 
p(TiGM|iHi)Pprior(iHi)/p(TiGM). We explicitly computed 
the term '"'j'^'^^^^l^^''' = In 10 x Tigm x piTiGM\xm) in 
this paper. Examples were shown in Figure 2 (Shi ~ 
0.91) and Figure 4 {xm ~ 0.80). The term Pprior(SHi) 
is the prior probability distribution for xm, which we 
assumed to be fiat. This expression only gives us the 
conditional PDF for ^hi for a given value of Tigm. 
To get a marginalized PDF, we compute p(a;Hi|TiGM > 
0.1) = J^^° p{xm\TiGM)p{TiGM)dTiGM, which simplifies to 
p{xm\TiGM > 0.1) = Ji^'° p{TiGM\xm)dTiGM (for our as- 
sumption that Pprior(a:Hi) ~ 1). 

If we require that Tigm > 0.1, then the blue dotted line 
shows that for the model without winds, a low neutral frac- 
tion is preferred. However, the statistical significance of this 
statement is weak: for this model xm < 0.70 at ^ 95% 
CL. Once winds are included, this constraint becomes even 
weaker, and it is not even possible to 'strongly' (i.e. > 95% 
CL) rule oulQ xui ^ 0.91. We caution against interpreting 
these numbers literally, as the overall shape of the PDF is 
completely dominated by the assumed prior on xm- This 
illustrates that it is difficult to rule out even large values of 
Xm with confidence-especially when winds are included in 
the modeling. 



4.2 Interpreting the Apparent Fast Drop in the 
'LAE Fraction' among LBGs 

IStark et al.l (|2010al ) recently found that the fraction of drop- 



^ Note that we switched to the shorter notation from probability 
theory, and denote probability density functions by p{a) = "^^^"^ , 
and conditional PDFs by p{a\b) = '^^^"^''^ 

^ If we require that Tigm > 0.3, then we can rule out xm ^ 0.80 
at ~ 95% CL for both wind models. Given existing uncertainties 
on both observed and intrinsic REW of the Lya line, requiring 
the IGM to transmit more than 30% is not well motivated. We 
therefore consider any constraints on xm from this requirement 
not relevant. 
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Figure 7. Wc show the CDF for REW. The solid line shows 
an exponential Lya REW-distribution that represents a simpli- 
fied representation of the z = 6 drop-out population (see text). 
The dashed line [dotted line] shows the z = 7 PDF under the 
assumption that the observed LAE fraction at z = 7 is different 
only because of the IGM, that the Universe was 51% [30%] neu- 
tral by volume (i.e. xm = 0.51 [Shi = 0.30]), and for the wind 
model with (Afuiifwind) = (10^'^cm ~^,25 km s~^). We also s how 
the upper limit that was derived bv lStark et al. I ll2010al , [201ll') for 
z = 7 drop-out galaxies. Explaining the observed rapid change 
in the LAE fraction among the drop-out population with reion- 
ization would require a fast evolution of the neutral fraction of 
hydrogen in the Universe. 



out galaxies with strong L yg emission decre ases strongly 
from z = 6 to z = 7 (also see lStark et al.|[201ll ). More speci- 
ficially, they found that the fraction of drop-out galaxies 
with a REW >75 A Lya line decreased by a factor of ~ 2 
be tween z = 6 and z = 7. A similar observation was made 
by lFontana et all (|201(]| '). who detected (weak; REW~ 13 A) 
Lya emission in only 1 out of 7 candidate z = 7 drop-out 
galaxies. These trends are not seen at z <^6. Such a sud- 
den, strong suppression of Lya flux from galaxies at z = 7 
could be a signature related to reionization, a lthough the 
available dat a still has large uncertainties. The IStark et al.l 
l|2010al . l201ll ) sample is consistent with no redshift evolution 
between 2 = 6 and z = 7 at the ~1 a l evel. And the s t atisti - 
cal significance of the drop observed bv lFontana et all (|2010l ) 
implicitly relies on the assumption that all seven drop-out 
galaxies are indeed at z ~ 7. Nevertheless, it is an interesting 
exercise to interprete these observations as is. 

In this work, we computed the Tigm-PDF at various 
ionization stages of the IGM during the EoR. Suppose that 
the observed REW-distribution at 2 = 6 is described by an 
exponential function (which provides a good fit to observed 
LAEs at 2 = 2-4, see Gronwall et al. 2007, Blanc et al. 2011) 
between REW= and REW= 300 A, i.e. Pz=6(REW) oc exp 
[-REW/REWc]. If we choose a scalelength of REWc = 50 
A, then the fraction of drop-out galaxies with REW ^75 A 
is f ^ 0.2, whi c h corre sponds to the median value observed 
by lStark et"al] (|2010al ) at 2 ~ 6. 

We now conservatively assume that the IGM at 2 = 6 
was 100% transparent to Lya photons emitted by galaxies, 
and that the observed REW-PDF at 2 = 7 is different only 
because of evolution of the ionization state of the IGM. Un- 
der this assumption, an observed REW at 2 = 7 requires 
an intrinsic (i.e. emitted) equivalent width of REW/TEgm, 
and we can compute the observed REW-PDF at 2 = 7 



as P.=7(REW) = AA/(3'd7iGMP(7!GM)Pz=6(REW/7iGM). 

Here, P(7i:GM) denotes the 7igm-PDF computed in this pa- 
per. The equation sums over all possible Tigm with the 
proper probabilities that a galaxy had a Lya line with 
REW/7igm md that the IGM transmitted a fraction Tigm- 
Finally, the factor Af normalizes Pz=7(REW). 

Figure [7] shows our model CDF of REW for 2 = 6 drop- 
out galaxies as the solid line, and for 2 = 7 drop out galax- 
ie^a;Hi ~ 0.51 (a;Hi ~ 0.30) as the dashed Ime {dotted Ime). 
We used the wind model with (A'hi, «wind) ~ (lO'^^cm"'^, 25 
km s~M- 

IStark et all (|2010al . boill ) put the upper Hmit on the 
fraction of drop-out galaxies with Lya REW^ 75 A at 
;^0.10 (also shown in Fig 7). The plots show that at in or- 
der to explain the observed evolution between 2 = 6 and 
2 = 7, we require a rapid evolution of the neutral fraction 
of hydrogen in the Universe (i.e. Aaini ~ 0.5 over A2 = ljf|. 
We stress that we assumed that Tigm = 1 at 2 = 6. If we 
had instead used a Tigm-PDF appropriate for a fully ionized 
medium at 2 = 6, then we would have needed the IGM at 
2 = 7 to be even more opaque, which required an even faster 
evolution in xhi- As we mentioned already in § [T] dust (as 
well as evolution in metallicity) would also require a faster 
evolution in Shi- For larger wind velocities and/or HI col- 
umn densities, we would again need a larger volume fraction 
of HI. Additionally, the data of lFontana et al.l (|2010l ) implies 
an even stronger evolution in the observed REW-PDF, and 
hence in the overall neutral fraction. 

Theoretically, the above-inferred rapid redshift evolu- 
tion in Xm is unrealistic even in models with no negative 
feedback on the source population (e.g., Barkana & Loeb 
2001; Fig. 9 in Mesinger et al. 2006; Fig. 1 in Lidz et al. 
2007). Furthermore, the sinks of ionizing photons (Lyman 
limit absorption systems) likely further slow the final stages 
of reionization (e.g., Furlanetto & Mesinger 2009; Alvarez 
& Abel 2010; Crociani et al. 2011), whose photoevapora- 
tion timescales could be much longer than A2 --^ 1 (Iliev et 
al. 2005). The inferred rapid evolution could mean that the 
current sample of 2 = 7 drop-out galaxies trace a region of 
the Universe that was more neutral than average. Alterna- 
tively, it could signal other physical effects: for example, a 
decreasing wind strength or covering factor towards higher 
redshifts can enhance the impact of the ionized IGM, and 

* Of course, wo computed the Tigm-PDF at 2 = 8.6, and now 
apply our models to 2 = 7 data. As a result, our model overes- 
timates the Gunn-Peterson optical depth. We therefore overesti- 
mated the damping wing opacity of the neutral IGM. If we had 
included the proper damping wing optical depth (i.e. at the cor- 
rect redshift), we would again have required an even larger xui 
to obtain the same total optical depth. Our current inferred rate 
of the evolution of Shi is therefore conservative. 
^ Our constraint on the redshift-evolution in xm depends some- 
what on the assumed range over which exponential function pro- 
vides a good fit to the data. For example, under the extreme 
assumption that there were no drop-out galaxies at z = 6 with 
REW^ 20 A, then we w ould need a s caleleri gth of REWc ~ 40 
A to be consistent with IStark et al.l l l2010al ). In this case, we 
would 'only' r e quire that Axhi ~ 0.3 over Az = 1. However 
IStanwav et al.l l l2007t) found that the fraction of drop-out galax- 
ies with weak emission (REW ^25 A) is consistent with the ob- 
served fraction at z = 3, which suggests that this 'extreme' model 
is unrealistic. 
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thus increase the rate at which the IGM opacity changes. 
Yet another interesting possibihty is that the observed rapid 
evolution in the Lya REW-PDF of the drop-out galaxy pop- 
ulation is entirely due to reionization, but that the Universe 
at 2 = 6 still contained a non-negligible volume fraction of 
neutral hydrogen (for a discussion on the current observa- 
tional constraints, see Mesinger 2010 and McGreer et al., 
2011). This is because the fraction of galaxies for which 
TiGM > 0.3 (as an example) evolves more rapidly at xm >0.5 
(i.e. df/dxm is largest when xm ~ 0.5; see Fig[5ll. For exam- 
ple, the fraction of galaxies with Tigm > 0.3 changes more 
rapidly between xui = 0.5 and xm = 0.7, then between 
xm = 0.0 and xm = 0.5 (see Fig [5]). This scenario how- 
ever would require an extended early epoch of reionization 
to be consistent with WMAP observations, perhaps driven 
by negative feedback on smaller-mass sources. Clearly it will 
be important to constrain the evolution of the LAE fraction 
between z = 6 and z — 7 with a larger sample of galaxies. 

Finally, we point out that the observed drop in the 
LAE fraction between z = 6 and z = 7 is consistent 
with observations of the LAE populations at z = 5.7 and 
z = 6.5. The Lya luminosity function of LAEs evolves sig- 
nificantly between z = 6.5 and z = 5.7 (IShirnasaku et al.l 
l2006l : iKashikawa et al.l l2006l : lOuchi et all |2010| ). However, 
the rest-frame UV luminosity function of these same galax- 
ies does not evolve between th ese redshifts (Kashikawa et 
al. 2006). iDiikstra et al.l (1200781 ) showed that these two ob- 
servations combined translate to a reduction in the num- 
ber of detected Lya photons from z = 6.5 by a factor of 
1.1 — 1. 8 (95% CL) relativ e to 2; = 5.7, with a median value of 
- 1.3 (|Ouchi et al.ll2010l '). Hu et al. (2010) recently obtained 
spectroscopic observations of narrowband selected LAEs at 
2 = 5.7 and 2 = 6.5 and found that the Lya REW of the 
average 2 — 5.7 [2 — 6.5] spectrum was REW~ 34 ± 2 A 
[REW~ 23 ± 3 A] . This corresponds to a reduction in the 
number of detected Lya photons from z = 6.5 by a fac- 
tor of ~ 1.3 relative to 2 = 5.7, consistent with the value 
inferred from the redshift evolution of the luminosity func- 
tions. This provides evidence for evolution in the REW-PDF 
th at is similar t o the t rends seen in the drop-out population 
bv lStark et"al] (|2010al ') and lFontana et al.1 (|2010l ). 



5 COMPARISON TO PREVIOUS WORK 

We found that Tigm ^50% for the majority of galaxies, even 
when the Universe is ~ 60% neutral by volume. In our mod- 
els, the IGM is more transparent than in previous studies of 
the visibility of LAEs during the EoR. Mesinger & Furlan- 
etto (2008b) assumed a constant opacity in the ionized IGM, 
and only considered velocity offsets in the Lya line due to 
the peculiar velocities of the dark matter halos (which are 
much smaller than the velocity offsets that can be imparted 
by winds). McQuinn et al. (2007) studied a very similar 
model, but did not include the halos' peculiar velocities. 
McQuinn et al. (2007) also investigated a 'wind model' in 
which the Lya line was redshifted by 400 km s~^, and found 
a considerable boost in the IGM transmission fraction for 
Xm ^0.7 (see their Figure 6). Iliev et al. (2008) studied the 
impact of the ionized IGM on the Lya line profile in more 
detail, but assumed Gaussian emission lines, which can re- 
duce the IGM transmission factor significantly compared to 



models that include outflows (see §|3TT]for a clear illustration 
of this effect). 

Perhaps surprisingly, our quoted values for Tigm are 
significantly higher than the values quoted in previous works 
even fo r a fully ionized IGM. For example, IDiikstra et al.l 
(|2007bD found that the ionized IGM could transmit as little 
as Ti:GM=0. 1-0.3 at 2 > 4.5 for a fully ionized IGM (also see 
Dayal et al. 2010). Laursen et al. (2011) recently c onstrained 
Tigm = 0.26t°Q ll at 2 = 5.7. IZheng et al.] l|2010al ) find that 
Tigm ~ 0.01 — 0.3, which in detail depends on luminositj(3- 

The main difference between the present and previous 
work are due to outflows, which causes most of the line 
flux to be systematically redshifted relative to the galaxy's 
systemic velocity (see §[3}. We acknowledge that the out- 
flow models were calibrated on the models of Verhamme 
et al. (2008) which assumed that the IGM at 2 = 2 - 5 
had no impact on the observed Lya line shapes. It is cer- 
tainly possible-especially towards higher redshifts-to repro- 
duce observed Lya line shapes with 'weaker' (i.e. lower ve- 
locity) outflows on ce the IGM is include d. This is illustrated 
by recent work of iLaursen et al.l (| 201lh. who modele d Lya 
RT in simulated, dusty galaxies. iLaursen et al.l (|201lh found 
that Lya photons emerge from their simulated galaxies with 
a broadened, double peaked profile (also see Barnes et al. 
2011). The IGM at 2 > 5 cuts off the blue peak, which then 
results in a typically observed redshifted, asymmetric Lya 
emission line (although at 2 = 3.5 the blue peak remains vis- 
ible for a significant fraction of sightlines through the IGM). 
In these models, galactic outflows have little impact on the 
Lya photons emerging from galaxies. 

While it is likely that Lya line shape is affected by 
the IGM to some extend (especially towards higher red- 
shifts), outflows are also expected to play at least an im- 
portant role. This is because outflows are ubiquitous in 
observed star forming galaxies, a nd the outflowing mate- 
rial has a large covering factor (|Steidel et al.1 bOlOl ). Im- 
portantly, there is evidence that outflows affect observed 
Lya line spectra: (i) asymmetric line shapes are present 
at low reds hift (2 ~ 0), when the IGM should not hav e 
an impact (|Mas-Hesse et all l2003l : iHeckman" et al.l I2OIII ): 
(ii) observations indicate that outflo ws promote the es- 
cape of Lya fro m a dusty medium (iKunth et al.l [l998l : 
Atck ct al.' '2008"); (m) 'Schaerer fc Verhammd (120081 ) and 
Dcssaugcs-Zavadsky ct al.. (20 10;) successfully reproduced 
the Lya spectrum of the Lyman Break galaxies cB58 and the 
'8 o'clock arc' respectively, with outflow mod e ls very sim i- 
lar to our own and those by IVerhamme et al.l l|2006l . |2008| ). 
but whose parameters were constrained by low-ionization 
metal absorption lines. In these models, the wind parame- 
ters inferred from the Lya line shape were consistent with 
those inferred from alternative observations. Steidel et al. 
(2010,2011) also constructed a simple model for the Lya 



Note that IZheng et al. ] lioioi) include scatte ring inside the 
virial ra dius in their calculat i ons, w hile for example lDiikstra et al.l 
ll2007bl) and ILaursen et al. I 1 I2OIII) do not. Gas within the virial 
radius can be significantly den ser and can have la rger peculiar 
velocities. This can explain tha tlZhe ng et a l.l ll2010al ) found lower 
values for Tigm- The fact that lZheng et a~ ri l2010al ) properly ac- 
count for the possibility that photons scatter back into the LOS 
(in close proximity of the source) can only boost their Tigm com- 
pared to that of other groups, and cannot explain this difference. 
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spectra observed from LBGs-as well as the Lya halos that 
are observed around LBGs-in which Lya photons scat- 
ter primarily through a large-scale galactic outflow, whose 
structure is constrained by low-ionization metal absorption 
lines. Points (i-iii) underline the probable physical connec- 
tion between the outflowing and scattering media. 

We have shown that our main conclusions remain valid 
for wind velocities in exces^^ of «wind ^25 km s~^ (and 
Nm ^ 10^*^ cm^^), which systematically redshifts the Lya 
emission line as it emerges from a galaxy. We consider these 
requirements to be reasonable, given direct observational 
constraints on o utflow velocities from low-ionization metal 
absorption lines jSteidel et al.|[201ol : iRakic et al.ll201(j ). and 
on HI column density (as in cB58 and the '8 o'clock arc', see 
e.g. Schaerer & Verhamme 2008 and Dessauges-Zavadsky et 
al. 2010),. 

After the submission of this work, Dayal & Ferrara 
(2011, hereafter DFll) submitted a paper in which they 
constrained xui < 0.2 using the same z — 8.6 galaxy. We 
argue that this upper limit is not robust for several reasons. 
Firstly, the model of DFll does not include outflows. In 
contrast, we have summarized the strong observational ev- 
idence that galactic outflows affect the Lya radiation field, 
implying that outflows must be included at some level in 
a realistic model. Secondly, there is an absence of large HII 
'bubbles' in the model of DFll: the HII bubble radii in their 
model are ;^3 — 4 Mpc when xui = 0.2, while simulations of 
reionization that include radiative transfer of ionizing pho- 
tons, show that the typical HII bubbles have radii that are a 
factor ~ 10 larger when xui ~ 0.3 (see e.g. Fig 2 of Zahn et 
al. 2010). Our large scale semi- numeric simulations of reion- 
ization properly capture the bubble size distribution, and 
the fact that the more massive galaxies preferentially popu- 
late these large HII bubbles. This absence of large bubbles in 
the models of DFll enhances the IGM opacity, which causes 
their constraints on a;Hi to be stronger than those obtained 
from our models where galactic winds are not included. 



6 CONCLUSIONS 

In this paper we have studied the visibility of the Lya 
emission line during reionization. We combine large scale 
semi-numerical simulations of cosmic reionization with 
empirically-calibrated models of galactic outflows. With 
these sophisticated tools, we compute the PDFs of the 
IGM transmission fraction, Tigm. We find that winds cause 
TiGM ^10% [50%], for the majority of galaxies, even when 
the Universe is 80% [60%] neutral by volume. This only 
requires wind speeds greater than ~ 25 km s^^, which are 
quite conservative judging by the observed Lya lines shapes 
at z < 5 (Verhamme et al. 2008, also see § [5]). Therefore, 
we conclude that the observed strong Lya emission from the 
reported z = 8.6 galaxy is consistent with a highly neutral 
IGM. 

We also show that evoking reionization to explain the 
observed drop in the 'LAE fraction' (see § 14. 2p of drop-out 



If wc lower A^Hi by an order of magnitude, then our main con- 
clusions are unaffected provided that ?^wind 

^200 km s-^ These 

requirements are also reasonable. 



galaxies between z = 6 and z = 7 (|Stark et al.ll2010al . l201lh . 
requires a very rapid evolution of Tigm, corresponding to 
Xui ~ — >■ 0.5 over Az = 1. Reionization models find such 
a rapid evolution unrealistic, which may indicate that either 
(j) the current sample of drop-out galaxies at z = 7 hap- 
pened to populate a region of our Universe that was more 
neutral than average, (ii) winds become weaker and/or have 
smaller covering factors towards higher redshifts, or (iii) that 
the Universe a.t z — 6 still contained a non- negligible volume 
fraction of neutral hydrogen. However, these conclusions are 
tentative as the available data still has large uncertainties. 

Regardless of these current observational uncertainties, 
our work underlines the point that Lya emission can be de- 
tected from galaxies in the earliest stages of reionization. 
This is a positive result for (narrowband) searches for high 
redshift Lya emitters such as the 'Emission-Line gal axies 
with VISTA Survey' (ELVIS) (e.g. iNilsson et al.ll2007l V On 
the other hand, if a neutral IGM is quite transparent to Lya 
photons, then a signature of reionization may be more dif- 
ficult to extract from observations of Lya emitting galax- 
ies. However, the redshift evolution of quantities such as 
(j) the 'LAE fraction'-or more generally the Lya restframe 
equivalent width PDF- among LBG s (|Stark et al.1 l2010al : 
iFontana et"alll20ld : [Stark et al.ll201ll 'l. and (ii) the UV and 
Lya luminosity functions of LAEs (Kashikawa et al. 2006), 
already provide interesting and useful constraints on mod- 
els of r eionization. Furtherm ore, the clustering signature of 
LAEs (|McOuinn et al.ll2007l : Mesinger & Furlanetto 2008b, 
though see Iliev et al. 2008) is also affected by reionization, 
and it has already been shown th at winds do not affect this 
prediction (|McOuinn et al.ll2007f) . 
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